Soft-Gluon-Pole Contribution in Single Transverse- Spin Asymmetries 

of Drell-Yan Processes 



J.P. and H.Z. Sang^ 

^ Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100190, China 
^ Center for High-Energy Physics, Peking University, Beijing 100080, China 
^ Institute of Modern Physics, School of Science, East China University of Science and Technology, 130 Meilong 

Road, Shanghai 200237, P.R. China 
^ ■ 

o 

(N 

P> i' Abstract 

' We use multi-parton states to examine the leading order collinear factorization of single transverse- 

spin asymmetries in Drell-Yan processes. Twist-3 operators are involved in the factorization. We find 
that the so-called soft-gluon-pole contribution in the factorization must exist in order to make the 
factorization correct. This contribution comes from the corresponding cross-section at one-loop, while 
(~| \ the hard-pole contribution in the factorization comes from the cross-section at tree-level. Although 

Oh' the two contributions come from results at different orders, their perturbative coefficient functions in 

the factorization are at the same order. This is in contrast to factorizations only involving twist-2 
D ■ operators. The soft-gluon-pole contribution found in this work is in agreement with that derived in 

i-C \ a different way. For the hard-pole contributions we find an extra contribution from an extra parton 

process contributing to the asymmetries. We also solve a part of discrepancy in evolutions of the twist- 
PsJ \ 3 operator. The method presented here for analyzing the factorization can be generalized to other 

^ ■ processes and can be easily used for studying factorizations at higher orders, because the involved 

\ calculations are of standard scattering amplitudes. 

(N 

(N 
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1. Introduction 

Single transverse-spin asymmetries(SSA) have been observed in various experiments, where an in- 
^ ' volved hadron is transversely polarized. A review about the phenomenologies of SSA can be found in [1]. 
%-{ . In general SSA can be generated if scattering amplitudes have nonzero absorptive parts and there are 
helicity-flip interactions. It has been shown in production of heavy quarks like top quarks there are size- 
able SSA[2l|3]. Because a top quark is heavy, the bound state effects can be neglected. The helicity-flip 
is due to the nonzero quark mass. Therefore, in the studies of [21 [3] one essentially deals with point-like 
particles and can use the perturbative theory in standard way. 

For SSA involving light hadrons, the origin of SSA is unclear because of bound-sate effects and the 
helicity-conservation of QCD with light quarks which can approximately be taken as massless. However, 
certain predictions can be made for SSA in cases where large momentum transfers are involved. In these 
cases, one can use the concept of factorization in QCD. SSA can be factorized in the form of a convolution 
with nonperturbative matrix elements of hadrons and perturbative coefficient functions. With this form 
of predictions one is then able to explore hadron structures with experiment. The coefficient functions 
can be calculated as a perturbative expansion. In this work we study the collinear factorization of SSA 
in Drell-Yan processes. We will show that the factorization of collinear divergences corresponding to 
nonperturbative effects is not made in an usual way as factorizations at leading twists. 
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The collinear factorization for describing SSA has been proposed in [11[5]. With the colhnear factor- 
izations SSA in various processes has been studied in [6l [71 [H [U [101 [HI [IS]- In such a factorization, the 
nonperturbative effects of the transversely polarized hadron are factorized into twist-3 matrix elements, 
or called ETQS matrix elements. Taking SSA in Drell-Yan processes as an example, SSA is factorized as 
a convolution of three parts: The first part is the standard parton distribution function of the unpolar- 
ized hadron defined with twist-2 operators. The second part consists of matrix elements of the polarized 
hadron defined with twist-3 operators. The third part consists of perturbative coefficient functions. The 
differential cross-section is determined by those twist-2- and twist-3 matrix elements of the initial hadrons 
and a forward hard scattering of partons from the initial hadrons. The perturbative coefficient functions 
describe the forward hard scattering. If the factorization can be proven, the coefficient functions can be 
calculated safely as a perturbative expansion and they are free from any soft divergence like collinear- 
and I.R. divergence. In this approach the effects of helicity-flip are parameterized with twist-3 matrix 
elements, while the absorptive part is generated in the hard scattering of partons. 

The above mentioned collinear factorization has been derived in a rather formal way by using diagram 
expansion at hadron level, in which one divides diagrams into three parts. Two of them are related to 
the initial hadrons respectively, the remaining one is related to the parton scattering. By expanding 
the two parts related to hadrons according to twist of operators and the part for the parton scattering 
with large momentum transfers, one obtains the factorized form of SSA and the perturbative coefficient 
function at leading order of as- The forward hard scattering is participated by three partons from the 
polarized hadron, i.e., two quarks with one gluon and two partons from the unpolarized hadron, e.g., 
q+{q-\-g)^j*-\-X^q + q with the antiquark from the unpolarized hadron and other partons from the 
polarized one. It seems difficult with this method to derive the perturbative coefficient function at higher 
orders and to prove the factorization. It is interesting to note that the contributions in the factorization 
consist of two parts. One part of the contributions is with the gluon carrying a nonzero momentum, 
called as hard-pole contributions, while another part of the contributions is with the gluon carrying 
zero momentum, called as soft-gluon-pole contributions. The two parts are associated with perturbative 
coefficient functions starting at the same leading order of a^. 

It should be noted that QCD factorizations, if they are proven, are general properties of QCD. These 
factorizations hold not only with hadron states but also hold when one replaces the hadron states with 
parton states. This is in the sense that the perturbatively calculable parts in factorizations do not depend 
on hadrons and are completely determined by hard scattering of partons. The procedure for this in the 
case of SSA is the following: With partonic states one can calculate the differential cross-section related to 
SSA and those twist-2 and twist-3 matrix elements with perturbative theory. In general they will contain 
soft divergences which usually appear beyond leading order. By writing the differential cross-section as 
a convolution of these matrix elements with a perturbative coefficient function, one can determine the 
function. Beyond the leading order of Og, one may be able to show that the function is free from soft 
divergences. If it is true, then the factorization for SSA is proven. This procedure also provides a way to 
determine the higher-order corrections to the perturbative coefficient function. 

In our previous works [TS j 114 1 115j '^^ have made such an attempt to derive the factorization by replacing 
hadrons with partons. To have helicity-flip with massless partons we have constructed a multi-parton 
state to replace the transversely polarized hadron in [15]. But with our partonic results at leading order 
of Os we can only find the hard-pole contributions. This is apparently in contradiction with the early 
results. 

In factorizations only with leading twist-2 operators, it is interesting to note the following fact: For a 
differential cross-section the factorizations at leading order is completely determined by partonic results 
at tree-level, i.e., the perturbative coefficient functions at leading order are determined only with the 
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differential cross-section and twist-2 matrix elements calculated at tree-level with parton states. This has 
the implication for one-loop results. If the factorization is right or proven, then the collinearly divergent 
part of the differential cross-section with parton states at one-loop is completely determined by the 
convolutions of the leading-order perturbative coefficient functions with the one-loop matrix elements of 
twist-2 operators. This patten about the collinearly divergent part of the differential cross-section can be 
iterated beyond one-loop. Assuming this is also the case for the factorization involving twist-3 operators, 
one then expects that the collinearly divergent part of the differential cross-section related to SSA at 
one-loop is completely determined by the convolution of the collinearly divergent parts of twist-2 and 
twist-3 matrices at one-loop with the perturbative coefficient functions at leading order. However, the 
assumptions may not be correct. Factorizations involving twist-3 operators can be different than those 
only with twist-2 operators. 

In order to clarify this issue we go beyond the leading order in this work. We find that certain con- 
tributions at one-loop, which contain collinear divergences, can not be factorized in the way as expected 
in the above. To make the factorization with the partonic states correct, one has to introduce additional 
contributions in the derived factorization which only contains hard-pole contributions. These additional 
contributions are just the soft-gluon contributions. This is an interesting fact because a part of the 
leading-order perturbative coefficient function is determined by quantities at non-leading order. It will 
be important for calculating higher-order corrections and proving the factorization, following the above 
outlined procedure. In this work we restrict ourself to the case where the forward hard scattering of 
those partons consisting of two antiquark from the unpolarized hadron and two quarks with one gluon 
from the polarized hadron. Beside the above mentioned contributions we also find an extra hard-pole 
contribution corresponding to the forward hard scattering where only the gluon is in the initial- or final 
state. By calculating the twist-3 matrix element corresponding to this hard scattering, we can solve a 
part of discrepancies in the evolution of the matrix element studied in [16^ \T7\ [TS] . 

In this work we study SSA in Drell-Yan processes in the kinematic limit of the small transverse 
momentum of the observed lepton pair. In this kinematic region there exists another factorization for 
SSA, called Transverse-Momentum-Dependent(TMD) factorization, similarly to the TMD factorization 
for unpolarized cases studied in |2H [22l [23l [2^ I25j . In the polarized case the nonperturbative effects of 
the polarized hadron are factorized into Sivers function[19] which contains both helicity-flip- and T-odd 
effects. The properties of Sivers function and SSA with it have been studied extensively [26l [271 129 [Ml 
[30l [311 [32l [33l [34l [35] . It should be noted that in the kinematic region of the small transverse momentum 
limit two factorizations apply, if the transverse momentum is much larger than the QCD scale Kqcd- 
It has been shown that the two factorizations are equivalent in the region [TOl [TTl I12j . Again, the TMD 
factorization here is derived in the formal way by using the mentioned diagram expansion at hadron 
level. In [131 [T4l [T5] we have examined the TMD factorization of SSA with parton states and found an 
agreement with existing results. In this work we will only focus on the collinear factorization. 

Our work is organized as the following: In Sect. 2 we give our notations for Drell-Yan processes and 
the factorization of SSA. In Sect. 3 we introduce our multi-parton states and give some relevant results 
of twist-3 matrix elements. In Sect. 4 we calculate SSA with our multi-parton state from certain classes 
of one-loop contributions which are collinearly divergent. We then show that these contributions can not 
be factorized as the usual way discussed in the above. These contributions in fact have to be identified 
as the mentioned soft-gluon-pole contributions in order to have finite corrections at higher orders of 
perturbative coefficient functions. In Sect. 5 we study the extra contribution which should be added 
to the factorization formula. There we also show that a part of the discrepancy in the evolution of the 
twist-3 matrix element derived in [16l[T71 is solved. Sect. 6 is our summary and outlook. 
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2. CoUinear Factorization of SSA in Drell-Yan Processes 

We will use the light-cone coordinate system, in which a vector is expressed as = {a'^ , a~ , a±) = 
((a° + a^)/\/2, (a° - a^)/\/2, a\ a^) and a^^ = (a^)^ + (a^)^. Other notations are: 

with the light-cone vectors I and n defined as Z'^ = (1,0,0,0) and n'^ = (0,1,0,0), respectively. We 
consider the Drell-Yan process: 

hA{PA, s) + hB{PB) -^7*{q)+X^r+i++ X, (2) 

where Ha is a spin- 1/2 hadron with the spin- vector s. We take a light-cone coordinate system in which 
the momenta and the spin are : 

Pa,b = {PXb^Pa,bM, = (0,0, sx). (3) 

Ha moves in the ^-direction, i.e., is the large component, hs moves in the — z-direction with P^ 
as the large component. The spin of hs is averaged. The invariant mass of the observed lepton pair 
is = (f'. The relevant hadronic tensor is defined as a matrix element of the forward scattering 
hA + hs^l* + X ^hA + Hb: 

W^'' = Y.j ^e^^-"(/iA(PA,s±),^B(PB)|g(0)7'^g(0)|X)(X|g(ar)7'^g(x)|/iB(PB),/i^(PA,sx)), (4) 
and the differential cross-section is determined by the hadronic tensor as: 



We are interested in the kinematical region where q\ <^ Q^. The hadronic tensor at leading twist 
accuracy has the structure: 

+q±a {eTQl + ^TqI) Q± ■ s±wP + ■■■ (6) 

In the above, we only give the structures symmetric in /jlu. W^\i = 1, 2, 3) represent T-odd effect related 
to the spin, wjj'"^^ are responsible for unpolarized cross-sections, wfj}^ contributes to SSA in the region 
q^ ^ which we will study. We introduce q~^ = xP^ and q~ = vPb- All structure functions depend 
on X, y and q"^. 

In the limit q± ^ only the structure function Wj!'^ gives the leading spin-dependent contribution 
to the differential cross-section and hence to SSA. The factorization of the structure function is the main 
subject to be studied in this work. In the coUinear factorization W^'^ is factorized with standard parton 
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distributions of hadron hs and twist-3 matrix elements of hadron Ha- There are two relevant twist-3 
matrix elements. They are defined as: 



TF{xi,X2)e'l s± 



2 



dyidy 



47r 



^-iy2(x2-X\)P+ -iyiXiP+ 



{{P.S^\i>{yinh^G+^\y2n)^mP,s^) - {s^ ^ -s^)] , 



dyidy2 
4tt 



'W2 



{x2—xi)P'^ —iyiXiP^ 



{{P, sx|V^(yin)7+75G+^(y2n)V(0)|P, s±) 



(7) 



In the above we have suppressed the gauge links along direction n between operators. These gauge links 
make the definition gauge invariant. The general properties of these twist-3 matrix elements have been 
discussed in [HI?]. One can show Ta,f{x,x) = 0. In general Tf{x,x) is not zero. This corresponds to 
the so-called soft-gluon-pole contributions because the gluon field in Tp{xi,X2) with xi = X2 = x carries 
zero momentum entering the hard scattering. 

For the case that SSA is generated through the scattering where an antiquark q from the unpolarized 
hadron hs and a quark q or a quark with a gluon from the polarized hadron h^, i.e, the forward parton 
scattering q + q + g^'y*+X^q + qoT the reversed, the structure function in the limit q±/Q <^ 1 can 
be factorized in the formllOl: 



Wi}\x,y,q^] 
A.h{x,yi] 
■As{x,yi] 



a. 



(27rg2)2 



dyi 
yi 



1 dy2 



TF{x,yi)- 



5(2/2) 
1-6)4 



Ahix, yi) + Asix, yi) + Asq{x, yi) 
+ TA,F{x,yi) 



1 



yi — (1 + 6)'^(1 - 6) +7>(yi,2/i) 



dyi 



- 3gf - 1 
, (1-6)+ 



'5(1-6) 



.M1±^<5(1 _ + 25(1 - 6)5(1 - ^2) In ^ 

(1-6)+ Q\ 



Asq{x,yi) 

6 



7>(yi,yi)5(i-6) 
x/yi, 6 = 2//y2- 



1 + 



6-1 

iV2 



1 + 6^ 
(1 - 6)+ 



25(1 -6) In ^ 



(8) 



In the above ^(2/2) is the antiquark distribution function of hs- The contributions to Wj^'' can be divided 
into three parts. The part with Ah consists of the hard-pole contributions. In Ah, the first term with 
Tf(xi,X2) has been first derived in [10] which is also confirmed in [15], while the second term with Ta,f 
has been derived in [15] . The part with As and that with Agq are of the soft-gluon-pole contributions, 
because they are related to Tf{x,x). The soft-gluon-pole contributions in Asq only appear in the limit 
in the limit Q ^ q±. 



3. Multi-parton State and Tvi^ist-S Matrix elements 

In [131 [H] we have studied SSA with single-parton states, where the helicity flip is caused be a finite 
quark mass m. In fact one can study SSA at parton level by taking massless limit m = 0. For this one 
needs to take the effect of helicity flip as that of a correlation between the spin of quarks and the spin of 



5 



gluons. For this purpose we consider the state or the system with the total hehcity A: 



with pi + k = p. In the first term 
specify the momentum as: 



n[X]) = \q{p,X,)[X])+cMPu^)9{k,Xgm), (9) 
A. For the qg-state, the total helicity is the sum Xg + Xg. We 



xqp'' 



p^ = 0, 0, 0), p'i = xop^, k^' = {l- xo)p^ 
The g-state and qg-state carries the same color index ic as given: 

\q{p,Xg)) = bl{p,Xgm, \q{puXg)g{k,Xg)) = T^.^blipi, XMi^, Xg 



(10) 



:iii 



where hi is the quark creation operator with i as the color index, ajj is the gluon creation operator with 
a as the color index, ci is taken real number. 

From standard text book we know that the transverse-spin dependent part of a matrix element, 
like the twist-3 matrix elements or the transverse-spin dependent part of W'^'^, corresponds to the off- 
diagonal part of the matrix element in helicity space. Because of helicity conservation in QCD with 
massless quarks, the twist-3 matrix elements are always zero if we replace the hadron in the matrix 
elements with a single quark. However, if one replaces the hadron with the above multi-parton state, 
the twist-3 matrix elements receive nonzero contributions from the interference between the single quark 
state and the state consisting of a quark and a gluon. In the interference, the quark always has the same 
helicity, while the helicity change is due to the helicity of the gluon. The structure function W^'^ also 
receives nonzero contributions from the interference. 

By replacing the hadron with the multi-parton state one can calculate those twist-3 matrix elements 
and the structure functions perturbatively for the purpose of factorization. At tree-level, it is straight- 
forward to obtain the twist-3 matrix elements as: 



Tf{xi,X2) 

7a,f(xi,X2) 



CigsTT\ 



{NI - 1)(X2 - Xl) [5{l - Xl)5{x2 - xo) - (5(1 - X2)5{xi - Xq)] , 

{N^ - 1)(X2 - Xl) [(5(1 - xi)(5(x2 - xo) + <5(1 - X2)<5(xi - xq)] . (12) 



2 ' ' 
/x"o 



These functions are proportional to ci indicating that they are from the mentioned interference. For 
simplicity we will set ci = 1 in the following sections without confusion. 
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Figure 1: A diagram for the twist-3 matrix element Tf{xi,X2) which gives a nonzero contribution with 
Xl = X2 = X. Black dots denote insertion of operators in the definition of rp(xi,X2). The broken line is 
a cut. 
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It is noted that at tree-level the function Tp{x,x) becomes zero. This is the reason why we can not 
re-produce in [15] the soft-gluon-pole contributions with tree-level results of Tp and Wj.^^ However, 
the function becomes nonzero at one-loop. To show this, we examine a particular contribution from a 
one-loop diagram given in Fig.l in Feynman gauge. The contributions contain a U.V. divergence and a 
coUinear divergence. Regularizing these divergences and subtracting the U.V. one we have: 



Tf{xi,X2] 



Fig.l 



Ta,f{xi,X2) 



Fig.l 



1 — Xi f Xq — Xi 




2x0 



-5{xi - xo) 



1 - Xo \ 1 - X2 



2x0 



(13) 



where the pole in ec = 4 — d represents the collinear divergence and fic is the scale associated with it. 
The scale // is associated with the subtracted U.V. divergence. The collinear divergence appears because 
the gluon going through the cut can be collinear to the incoming gluon and the outgoing quark in Fig.l. 
In the above we also give the contribution to Ta,f from Fig.l. From this result we see that Tp is nonzero 
at xi = X2- After examining all one- loop diagrams in Feynman gauge, we find that only the diagram in 
Fig.l gives nonzero contribution at xi = X2, i.e., 



7>(x,x) 



N,{N^ - l)xo^ 



+ 7 - In 



(14) 



This result is in agreement with our previous calculation in the light-cone gauge n ■ G 
function Ta,p is always zero at xi = X2 = x. 



0[T5]. The 
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Figure 2: The diagram for the tree- level contribution to SSA in the limit q± Q. The black dots 
represent the insertion of electro-magnetic current operators. The broken line is the cut. The short bar 
cutting the quark propagator means to take the absorptive part of the propagator. 

To calculate the hadronic tensor, we replace the polarized hadron with the parton state \n) 
of Eq.([9]), and the unpolarized hadron hs with an antiquark with the momentum = (0,p~,0,0). 
At tree-level, also receives contributions from the interference, i.e., from the forward scattering 

q + q + g^^*+g^q + q and q + q^j*+g^q + q + g. In the limit q± <^ Q, only one diagram 
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given in Fig. 2 gives the contribution. In the diagram the short bar means to take the absorptive part of 
the cut propagator with the momentum kq, i.e., 

= n-kq [-i7r6{k''q)) . (15) 

In fact, the short bar here represents a physical cut of the amplitude of the left part in the diagram. It has 
been shown in [15] with the tree-level results of the twist-3 matrix elements given in Eq.(12) that the tree- 
level result of only produces the Ah-tevm in Eq.dH]). If one expects that the factorization involving 
twist-3 operators here happens in the same way as in the factorization only with twist-2 operators, one 
will conclude that at leading order of as of the perturbative coefficient function is predicted only 
with the hard-pole contributions of Eq.([8]). This is obvious in contradiction with the results in Eq.Q. 
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4. Soft-Gluon-Pole Contributions 

At tree-level, all momenta carried by gluon lines in Fig. 2 are fixed and can not be zero, e.g., the 
momentum ki of the gluon crossing the cut is fixed by the total momentum conservation. Therefore, we 
can not identify any gluon line in Fig. 2 corresponding to the gluon field with zero momentum in Tf(x, x) 
in Eq.(7). Now we consider the case in which there is an extra gluon exchanged and crossing the cut as 
those diagrams given in Fig. 3. In this case the momentum ki has to be integrated. In the integration 
the collinear region, in which the gluon with ki is collinear to the incoming gluon and outgoing quark, is 
included. This will result in a collinear divergence. Taking Fig. 3a as such an example, there is an extra 
gluon exchanged in comparison with Fig. 2. The momentum ki hence will be integrated. In the collinear 
region of ki, one can realize that the part of Fig. 3a including the three gluon vertex and the vertex 
absorbing the gluon with ki is essentially given by Fig.l. This indicates that the collinear ly divergent 
part from the collinear region may be factorized with Tp{xi, X2) or Ta^f given in Fig.l in the form of 
a convolution of Tp or Ta,f with a perturbative coefficient function. If the function is not the same as 
those in Ah of Eq.dSj) determined at leading order, then one has to add extra terms beside the term with 
Ah in Eq.dS]) to make sure that the factorization is correct at one-loop level. It is interesting to note 
that by taking the absorptive part of the quark propagator in the left part of Fig. 3a, one finds that in 
the collinear region of ki the momentum of the gluon exchanged between the initial gluon and initial 
g is a soft gluon. More precisely, the exchanged gluon is a Glauber gluon with the momentum patten 
k'^ ~ (Aq, Ag, Aq, Aq) with Aq ^ 1. This indicates that the possible extra terms may be soft-gluon-pole 
contributions factorized with Tf{x,x) from Fig.l. In this section we show that this is indeed the case. 

As discussed in the above, it is easy to find those diagrams at one-loop which can give contributions 
to the soft-gluon-pole contributions. The diagrams are those where the incoming gluon emits a collinear 
gluon and the collinear gluon is absorbed by an outgoing collinear quark. Beside these gluons, there is an 
extra exchanged gluon crossing the cut responsible for the finite q±. A class of those diagrams is given in 
Fig. 3. For our purpose we consider those contributions from the collinear region where the gluon crossing 
the cut and emitted by the initial gluon is collinear to the initial gluon and the outgoing quark. 

We take Fig. 3a as an example to show how we obtain the collinearly divergent part. The contribution 
from Fig. 3a to the hadronic tensor can be written with standard Feynman rule as: 



jd_fa^a3a f d^d^^ 



2Ar/™y ^^5^(p + p-fci-A:2-g)(27r5(fcf))(2vr5(fc2'))Hvr5((A;3+p)^)) 
l;(p)7^^7 • (^3 + Ph^l • iPi - k2)Y'T^'T''^n{pi,\q) 
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Figure 3: The diagrams for the ampHtude q + q + G — )• 7* + G + G — )• g + which gives a part of 
contributions to SSA. 



■u{p, Xg)T"-^T''^jp^'j ■ {p - /ci)7p27 -{p-ki- k2)j''v{p) 
(e • {-ki + ks)g%, + {-ks - ky^e^, + {k, + k)^,eP' 



{{p - ki - k2Y - ie){{p - kiY - ie){{pi - k2Y + ie){kl + ie) ' 



(16) 



the color and spin of the initial antiquark q is averaged and it gives the factor (2A'^c)^^- The initial 
quark has the helicity \q, the initial gluon has A^. The absorptive part in the scattering amplitude is 
generated by the cut cutting the quark propagator. This gives the (5-function 5{{k2, +pf') with k^, being 
the momentum carried by the gluon exchanged between the initial gluon and the initial antiquark in the 
left part of Fig. 3a. We will consider the collinear region where the momentum ki of the gluon emitted 
by the initial gluon is collinear. In the collinear region the momentum ki scales as: 

fcf ~(1,A2,Ao,Ao) (17) 

with Ao ^ 1. The on-shell condition from the quark propagator fixes ki in the collinear region as 

^ -m5{{k^+pf) ^kt ^k-^ -^ + 0{\l). (18) 



{k^+py + ie ^ J i 2p- 

This constraint leads to that the gluon with k^ is a Glauber gluon. It is soft and may be represented by 
the gluon field in Tf{x,x). 

The evaluation of the contribution containing the collinear divergence from the collinear region given 
by Eq. ()17p is rather straightforward. One first uses these (5- functions to perform the integration over k2, 
k^ and kf. The remaining integration is that over ki±. The integrand is then besides some trivial factors 
a product of those terms in [• • •] in Eq. (fT6]) . the denominators of propagators and the 5-function (^(/cl). 
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Now one can expand the integrand in Aq- The leading order is at Aq ^ which does not give the cohinear 
divergence. The next-to-leading order is at Aq , which give the collinear divergence after the integration 
over ki±. Contributions from higher orders are finite. In the expansion we notice that the 5- function 
(^(/cl) also depends on ki and needs to be expanded. The expansion will give a contribution proportional 
to the derivative of the 5- function. This contribution may correspond to those terms in Eq.([8]) with the 
derivative of Tp(x,x). 

After the integration over ki± one can take the limit q± <^ Q. To derive the limit we will use the 
following in the limit of q± — )> 0: 



s5(s(l - y){xo -x) - qi) 
s6'{s{l - y){xo -x) - q\) 



'5(1 - y) , - x) q]_ 
[xo - x)+ (1 - y)+ 



11 



du 



(19) 



with 



2p^p . The calculation can be simplified by the following: We need only to calculate the 



off-diagonal part of the matrix element in helicity space. In this part one always has AgA^ 



-1. We win 



set XqXg 



— 1 in our calculation for simplicity. We have the collinearly divergent part of the hadronic 
tensor in the limit q± <^ Q : 



1 iiuA\) ■ 1± 



3a 



(4vr) 



Nr. 



-91 



il 



2 \2 



Xq 



+ ••• 



(20) 



where • • • stand for the following contributions: the contributions at non-leading order in the expansion 
in q±/Q, the contributions which do not contain the collinear divergence and the contributions of tensor 
structures other than g'^ . These contributions are irrelevant for our purpose. By adding the complex 
conjugated contribution with different parton helicities as a part of the interference, one can then obtain 
the off-diagonal part of the hadronic tensor in the helicity space. From the off-diagonal part we can 
extract the structure function: 



Wr 



(1) 



3a 




y) 



+ ••• 



(21) 



again, • • • denote those irrelevant contributions which do not contain the collinear divergence or are not 
at the leading order in the expansion in q±/Q. In the limit the leading order of is of qj^. In 

the following we will only give the collinearly divergent contributions in the limit of gj_ — )■ explicitly. 
Performing similar calculations we have the contributions from other diagrams in Fig. 3: 



Wr 



(1) 



Wr 



(1) 



Wr 



(1) 



36 



3c 



3d 



0{ql'), 



(47r)2 N, 



^2x^ f xoyS{xo - x) ^ x5{l - y) 

{xq - x)+ 



/ ^2 \ 



+ 



+xq5{xo - x)5{l - y) I 1 - In ^ ) ^ + 



9sal - 
(47r)2 N, 



y/2x^ \ XQy'^5{xQ -x) ^ x^(5(l - y) 
{q\f\ (l-y)+ xo{xo-x)+ 



+xo5(xo - x)6{l - y) 1^1 - In 1^ j I + • • • , 



(22) 



10 



(«) 




Figure 4: The diagrams for the amplitude q + q + G 
contributions to SSA. 



7* + G + G— s-g + Q', which gives a part of 



Beside those diagrams in Fig. 3, there is another class of diagrams which give the wanted contributions. 
The calculations of these diagrams are slightly different than those of Fig. 3. We illustrate this by taking 
Fig. 4a as an example. The contribution from Fig. 4a can be written as: 



4a 



9 s 
2N, 



(27r)4 (27r)4 



5\p + pi -k2- q){2^5{kl)){m5{{h + pf) 



■v{p)j^j ■ {pi - k2)jp2T"-^T"-u{pi) (e • {-ki - k?,)gp^^.^ + (^3 - k)p^ep^ + {ki + k)^^ep^) 



{{p -ki- /c2)2 - ie){{p - kiY - ie){{pi - k2Y + i£){kl - i£){kl - ie) ' 



(23) 



Again, the initial quark has the helicity A^, the gluon has \g. The cut of the quark propagator gives 
the (5- function <5((A;3 with ^3 = /ci — /c as the momentum of the gluon emitted by the antiquark 

in the right part of Fig. 4a. ki is the momentum of the gluon emitted by the quark in the right part of 
Fig. 4a. Unlike the gluon with ki in Fig. 3a, where it is on-shell, the gluon with ki in Fig. 4a is off-shell 
in general. This difference results in that the integration over fcj" looks nontrivial at the first step, while 
the integration over /cf in Fig. 3a can be simply done with the on-shell condition 5{k\). 

Now we consider the collinear region where ki is collinear to k and p. The scaling of its each 
component is given in Eq. (jl7p . Then from the (5-function 5((/c3 -|- p)^) kf is fixed as k^ = k'^ + ©(Ag) 
after the integration over k^ . The integration over k^ can be done with a contour in the complex k^- 
plan. With the fixed ki one can show from Eq. ()23p that there are poles from denominators of quark 
propagators only in the lower-half plan. These poles are corresponding physical cuts. One can use these 
poles by taking a contour in the lower-half plan to perform the integration. However, we notice that there 
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is only one pole in the upper-half plan. The pole is from the gluon propagator with the momentum ki. 
One can cquivalcntly use this pole by taking a contour in the upper-half plan to perform the integration. 
Therefore, the integration over A;f can be done easily by the replacement: 



%£ 



2Tri5{kf). 



(24) 



This also applies for other three diagrams in Fig. 4. It is interesting to note that the gluon with ki or 
with ks in Fig.4a are in corresponding to the gluon with ki or with ^3 in Fig.Sa, respectively. The gluon 
with A;3 in Fig.4a is also a Glauber gluon. The remaining calculations are similar to those of Fig.3. We 
have the following results of Wj^^ from Fig.4: 

^'2 Ar2 - 1 



Wr 



T |4a 



Wr 



T |46 



Wr 



T I4d 



(47r)2 N, 

9sO?s - 1 
(47r)2 Nc 

(47r)2 N, 
gsal iV2 _ 1 



(47r) 



Nr. 




x{xq - x)5{l - y), 
-y)6{xo -x), 
xoy5{xo - x), 



-\xd{l-y). 



(25) 



Summing the contributions from Fig.3 and Fig.4 together we obtain the coUinearly divergent contribution 



of W^^ , denoted as Wr 



(1). 

T,s- 



Wr 



(1) 



T,s 



(47r)2 
+<5(1- 



Nr 



C) 



(i-y) 




2<^(l- 0(^(1 -y) 
2(3- 3(2 _i 

(1 - 0+ 



+ 2C 




c 



X 

Xo' 



(26) 



If the factorization here takes the same patten as those only with twist-2 operators as discussed in 
the introduction and assuming that there is only the hard-pole contribution at leading order, one then 
expects that the above W^^l should obey: 



Wr 



(1) _ 



a. 



T,s 



(2vrc?2) 



>S{l-y) £ 



dyi 
vi 



:i-6)- 



-TF{x,yi] 



Fig.l 



+ TA,F{x,yi) 



(27) 



so that the coUinear divergences caused by the collinear gluon in Fig.3 and Fig.4 will not appear in the 
perturbative coefficient functions at one- loop. In the above we have already taken the tree-level result 
^(2/2) = (5(1 — 2/2)- It is easy to see that the above equation does not hold because the color factor of 
Wj^l does not match that of Tf(x, yi) and Ta,f{x, yi) from Fig.l. Therefore, the factorization at leading 
order must contain extra terms besides the hard-pole contributions. 

With the discussion at the beginning of this section, parts in each diagram in Fig.3 and Fig.4 can be 
identified with Fig.l. By deleting these parts, these one- loop diagrams reduce to those for the forward 
scattering q + q* +g* ^ +X ^ q + q*. In the kinematic region considered here, the off-shell quark can 
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be approximately taken as an on-shell quark. The virtual gluon is the mentioned Glauber gluon. It carries 
the momentum with the +-component approaching to zero. Therefore, W^l should be factorized with 
Tf(x,x) from Fig.l. With this observation, we realize that the factorization of this part is exactly given 
by the term with As of Eq.dS]). It is interesting to note that the term with the derivative of Tf(x,x) is 
also reproduced. Therefore, we can already conclude here that factorizations involving twist-3 operators 
are different than those factorizations only with twist-2 operators. The perturbative coefficient function 
in a factorization involving twist-3 operators can not be completely determined by tree-level results of 
the differential cross-section and twist-3 matrix elements. 

Besides the contributions studied in the above there are soft-gluon-pole contributions appearing in 
the limit q± — t- 0. There are two classes of diagrams given in Fig. 5 and Fig. 6. Their calculations are 
similar to those in Fig. 3 and Fig. 4. We only list their divergent contributions: 




Figure 5: The diagrams for soft-gluon-pole contributions appearing in the limit g_L — )• 



Wr 



T 1 5a 



N,{N^ - 1)^^ (-f ) V2^xo6ix - xo)y + 

2^S{x -Xq) 



(47r)2 


Nc \ 




(N^ - If 


(47r)2 






(iV2 - 1)2 


(47r)2 
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y 



(i-y)- 

-(5(x - a;o)(l - y), 
V 



-6{x - Xq) 



ql 



the contributions from Fig. 5b and Fig. 5c are non- leading in the limit of q±/Q <^ 1. 
The contributions from Fig. 6 are: 



Wr 



(1) 



g^al NcjN^ - 1) V2^ 

6a (4vr)2 2 (g2)2 

, <5(i-y) 



6{xo - x) 



xoy{2y - 1) 



{xq - a;)+ xq 



xo6{xo - x)6{l - y)ln 



qV 



(28) 
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Figure 6: The diagrams for soft-gluon-pole contributions appearing in the hmit g_L — )• 



(1) 



6b 



Wr 



(1) 



Wr 



T 
(1) 



(47r)2 2 
^ (5(1 - y) a:(23;o - x) 

(xq - x)+ Xo 
gsal - 1 



6{xo - x) 



(1-2/)- 



T 



6c 



6cZ 




xo6{xo -x)6{l -y)ln^2 
2;o2/(l - y)S{xo - x), 



xoySixo - x) 



(29) 



the contributions from Fig.6e and Fig.Gf are non-leading in the hmit of q±/Q <^ 1. In Fig. 6a the gluon 
with ki corresponds to the gluon with ki in Fig. 4a. We have observed that the integration over k^ here 
in Fig. 6a is different that in Fig. 4a. In the lower-half complex fcj^-plan there is only one pole from the 
quark propagator in the right part of Fig. 6a, and there are three poles from the three gluon propagators 
in the upper-half complex fc^-plan. One may take a contour in the lower-half complex /cj~-plan to perform 
the /cj~ -integration whose result is only from the pole of the quark propagator. By taken a contour in 
the upper-half complex k^-plan, this integration result can also be written a sum of contributions from 
the three poles of the gluon propagators. In the limit q±/Q <^ 1, the terms in Wj!'^ proportional to 
6{xq — x) only come from the contribution of the pole in the gluon propagator with ki, i.e., these terms 
can equivalently be calculated with the fc^-integration by taking Eq.(24) as for Fig. 4a. With this fact the 
gluon exchanged between the two gluon lines is a Glauber gluon. This indicates that these terms may be 
factorized with Tp(x, x) according to the experience from Fig. 4. The above discussed also applies for the 
remaining diagrams in Fig. 6. In Eq. (j29p the terms proportional to 6{1 — y) except the terms containing 
log's can be identified as hard-pole contributions at one-loop. 
The sum of the diagrams in Fig. 5 and Fig. 6 is then 



(1) 



Fig.b+FigA 



(4vr)2 (g 



2 \2 



-25(1 -y) In ^ 



+ 



Nc{Nl-l){5{xQ-x) 
x6{l - y) 



i + y' 



1 



y 



1 



Xq{xq 



(30) 
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With Tp(x,x) in Eq. ()14p it is clearly that the terms in [• • •] reproduce the term Agq in Eq-dS]). The 
last term in the above is not a soft-gluon-pole contribution and will become relevant if one studies the 
perturbative coefficient functions at the next-to-leading order. 

Before ending this section an interesting observation can be made. We observe that SSA calculated 
here at one-loop is generated through exchange of Glauber gluon at one-loop and it is divergent. This is 
in contrast to factorizations only with twist-2 operators for Drell-Yan processes. These factorizations are 
for those differential cross-sections which do not contain T-odd effects. In proving these factorizations, 
the existence of Glauber gluons brings up the most difficult obstacle[36l EH [38] • But it is able to show 
that the divergences caused by Glauber gluons are canceled in differential cross-sectons[36l EZl [38]. For 
the factorization studied here, such divergences are not canceled and they need to be factorized into the 
twist-3 matrix element with xi = X2- This will have some implications for the study of factorizations in 
the framework of soft collinear effective theories of QCD[39j. 



5. Additional Contributions 

In the previous sections we have used the multi-parton state in Eq.dOj) to replace the polarized hadron 
and determined the factorization form of Wij^^\ After the replacement, we have studied SSA essentially in 
the partonic forward scattering process q + [q + G) — )■ -)* + X q+q or the reversed scattering, where the 
helicity difference between the initial gi^-state and the final g-state is ±1. For a real hadron scattering, it 
is possible that instead of the above scattering one has the forward scattering q + {q + q) — >■ 'y*+X — )• q + g 
or the reversed, where the final gluon and initial qq come from the polarized hadron. If the total helicity 
of the qq state is zero, this forward scattering will also delivery an additional contribution to SSA besides 
those given in Eq.dH]), because the helicity difference in the scattering is also ±1. This has been realized 
in the study of the evolution of the twist-3 matrix element Tp(x. a;)|17j. 

The additional contribution can be factorized with the twist-3 matrix elements. The factorization 
can be studied with our multi-parton state in Eq.Q by adding a qqq state as a component: 

M^]) = k(p,Ag)[A]) + ci\q{pi, Xq)g{k, Xg)[X]) + C2\q{ki, Xi)q{k2, X2)q{k3, X3)[X]), 

fcf = Zip, for i = 1,2,3. (31) 

Replacing the transversely polarized hadron h^, one obtains SSA from the interference of the state qqq- 
state with the g^-state. The interference with the single quark state will not contribute to SSA because 
of the helicity conservation. By taking one quark from the qqq-state and g^-state as a spectator quark, 
one can have the mentioned forward scattering q + (q + q) ^ ^* + X ^ q + g oi the reversed scattering. 
Therefore, we need effectively to replace, e.g., in the twist-3 matrix elements the state with a qq 
state and the state (/i^l with a gluon, or in a reversal way. The total helicity state of the qq state should 
be zero. The qq state is in color octet in correspondence with the gluon. Taking the quark with ^2 as 
the spectator, one can simply work out those twist-3 matrix elements at tree-level: 

TF,c2ixi,X2) = -C2NTTgs{Nl - 1)xq{x2 - X\)^2zxZj, 

■ [6{xi + Z3)6{X2 - Zi) - 5{X2 + Z3)6{xi - Zi)] , 

TA,F,c2ixi,X2) = 0. (32) 

The factor A/" is a normalization factor of the spectator quark state with other trivial factor. The same 
factor will also appears later in and Tp{x,x). All relevant quantities calculated with the state 
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Eq. ()3ip are the sum of contributions studied in previous sections and those studied here. Because of this 
our study of these contributions can done separately. It should be noted that there are two possibilities 
to have one quark as a spectator for the interference with the state in Eq. (|3ip . For simplicity we only 
present the results of the above possibility in this section, because this will not affect the factorization, 
i.e., the perturbative coefficient functions. We denote quantities calculated in this way with an index C2- 
We have performed detailed calculations for these two possibilities and obtained the same results for the 
factorization form of W^'^ and for the evolution of Tp{x,x) which will be discussed later. 



P 




Figure 7: The tree-level diagram contributing to SSA. 



Similarly, when we do the replacement of for the hadronic tensor as done for the twist-3 matrix 
elements in the above, one can obtain nonzero SSA. At tree-level, in the limit q± <C Q, only one diagram 
gives nonzero contribution. Calculating the contribution in the similar way as shown in [15] and previous 
sections, we obtain the structure function: 



(1) 



C2M 



^''xo\/2zi 2:3^(1 - y)5{x - zi) ^ . 

X + Z3 



Comparing with the above Tp we can obtain a factorized form as: 



w!j^l^{x,y,qA_) 



a. 



C — t ^q{y2)TFix, X - yi)5il - 6)(1 - 26), 

Jx yi Jv y2 



(33) 



(34) 



it should be noted that the contribution of Tf{xi,X2) calculated in Sect. 3 will not be involved here 
because this contribution is zero for X2 < 0. It is clear that this part should be added to the factorized 
from in Eq.([8]), i.e., the term Ah should be modified as: 



Ah{x,yi) = 6(1-^2) 



Tpix.yi) ,} +T^,F{x,yi) + ^--^Tf{x,x - yi] 



iV2 



(35) 



(1-6)+ 

In the above we note that one argument of Tf{xi, X2) with X2 = x — yi'is negative, representing the fact 
that there is an antiquark from the polarized hadron. 

In [161 [IZj ths evolution of the twist-3 matrix elements have been studied. The results are different. 
The evolution of the non-singlet part of Tp(x,x) with x > is given with z = x/^ in [17] as: 



dTpix, X, fi) 



a. 



vr 



1 (ie 
J 



P,,{z)TF{i.i) + ^ 



Nj{l + z)Tf{x, - (1 + 2')7>(6 



1 - z 



+ 7a,f(2;,C) 



N,T{x, x) + ^J^ ^ [(1 - 2z)Tf{x, x-0- rA,F(x, x - 0] |. 



(36) 
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The discrepancies is the following: From [16] the evolution has only the terms in first line. The last term in 
the first line has a different sign than that of [17]. With our multi-parton state one can calculate Tp(xi, X2) 
to check the evolution. For this we calculate Tf(x,x) with the contribution from the interference with 
the qqq state in the same way showing in the above. 



Figure 8: The one-loop diagram to Tp^c-iix^x). 
For X > there is only one diagram at one-loop showing in Fig. 8. We obtain the following: 



Tf,c2[x,x,h) = C2N — XQ^/2yiyz6{x -yi) — ■ 

4A'c x + y^ 



2 \ 

7-Mn47r + In ^ 

ec J K 



■ (37) 



Adding the results from Sect. 3. we have: 



dTp{x, X, fi) 
din /J, 



as 
+ • 



1 

J 



Nr. 



1 + z 



T>(x,0 +TA,Fix,£,) 



2z)TFix,x-0 



(38) 



Comparing the above result with that in Eq. (j36p . we find an agreement except the terms with Tp(x,x) 
and the last term in Eq. (j36p . These terms are denoted with • • • in Eq. (j38p . They can not be obtained 
at one-loop with our results here. To verify them one has to study them at two-loop level. With our 
one-loop results at least a part of the discrepancy is solved. 

From the interference with the qqq state one also expects that there are soft-gluon-pole contributions 
in Eq. (|34p . These contributions are generated in a similar way as those studied in Sect. 4. Wj^^^ receives 
contributions at one-loop level from similar diagrams in Sect. 4., where one replaces the incoming gluon 
line going through the cut with an antiquark line and the outgoing quark with an outgoing gluon line. 
Comparing the topology of these diagrams with that of Fig.8 for Tp^c2{xi x), one can expect that the 

contributions to Wj^l^ from those diagrams can be factorized as the soft-gluon-pole contributions in 
Eq.([8|). We have checked these contributions. The results confirm the above expectation. 



6. Summary and Outlook 

Because of the helicity-conservation of QCD, SSA of a single quark state can not be generated in the 
involving forward scattering. We have used multi-parton states in order to have a nonzero SSA in the 
relevant partonic processes. Using the multi-parton states one can calculate SSA of Drell-Yan processes 
and relevant twist-3 matrix elements. With these partonic results we can examine and derive the collinear 
factorization. The collinear factorization for SSA has been derived in the literature based on a diagram 
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expansion of the relevant hadronic tensor. By using partonic results at tree-level, only the hard-pole 
contributions can be obtained. The soft-gluon-pole contributions can not be obtained. The reason for 
this has been discussed. 

In this work we have performed the study with multi-parton state at one-loop level in order to ex- 
amine and identify these soft-gluon-pole contributions in Drell-Yan processes. If we assume that the 
factorization is correct and there are only hard-pole contributions at leading order of perturbative coef- 
ficient functions, we find that a class of one-loop contributions to SSA, which are collinearly divergent, 
can not be factorized at one-loop order. To correctly factorize collinear divergences appearing in these 
contributions, one has to add extra terms in the factorization derived with tree-level partonic results. 
Interestingly these extra terms are just those soft-gluon-pole contributions. Therefore, with our multi- 
parton states we can re-derive the soft-gluon-pole contributions which have been derived with the diagram 
expansion in[10j. 

It is interesting to note that the hard-pole- and soft-gluon-pole contributions in SSA, i.e., in the 
structure function are at different order of Os- But, the perturbative coefficient functions are at 

the same order. The perturbative coefficient functions associated with the soft-gluon-pole contributions, 
although derived from SSA at one-loop, are at the same order of those associated with the hard-pole 
contributions, which are derived from SSA at tree-level. This is in contrast with factorizations for differ- 
ential cross-sections, where only leading twist-2 operators are involved. In these twits-2 factorizations, 
perturbative coefficient functions at leading order are completely determined by differential cross-sections 
at tree-level and tree-level matrix elements of involved twist-2 operators. Therefore, our study here also 
shows an unusual feature of factorizations involving twist-3 operators. 

By taking multi-parton state we also find a new contribution to W^^^^ The new contribution comes 
from the parton process where a qq pair with the total helicity A = is transmitted into a gluon. This new 
contribution can be factorized with the twist-3 matrix element. One important twist-3 matrix element is 
calculated at one-loop in this work. From the result one can derive the evolution of the matrix element. 
In this work we can solve a part of discrepancy between evolutions derived in [16\ I17j . To solve the 
remaining parts one has to calculate the twist-3 matrix elements with the multi-parton state at two-loop. 

In this work we have restrict ourself to the certain relevant partonic processes with one antiquark from 
the unpolarized hadron. Having succeeded to reproduce the soft-pole-gluon pole contributions with multi- 
parton states, one can start to analyze with the method presented in this work other relevant partonic 
processes. E.g., the processes involving one gluon from the unpolarized hadron. In such processes, it is 
possible to have soft-quark-pole contributions represented by Tf{0,x). One can also use our method to 
analyze another type of SSA appearing in the case when the momentum of each lepton is measured. In 
this case the factorization of SSA is so far still different derived from different works [30]. In principle, 
our method is not restricted to Drell-Yan processes. It can be used in any case where SSA appears. We 
believe that it has advantages to use our method for analyzing factorizations of SSA and for calculating 
higher order corrections, because the involved calculations are of standard scattering amplitudes. 
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